Differential equations for deformed Laguerre polynomials 
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Abstract. The distribution function for the first eigenvalue spacing in the 
Laguerre unitary ensemble of finite size may be expressed in terms of a solu- 
tion of the fifth Painleve transcendent. The generating function of a certain 
discontinuous linear statistic of the Laguerre unitary ensemble can similarly be 
expressed in terms of a solution of the fifth Painleve equation. The method- 
ology used to derive these results rely on two theories regarding differential 
equations for orthogonal polynomial systems, one involving isomonodromic 
deformations and the other ladder operators. We compare the two theories by 
showing how either can be used to obtain a characterization of a more general 
Laguerre unitary ensemble average in terms of the Hamiltonian system for 
Painleve V. 



1. Introduction 

1.1. Objective. Spacing probabilities and moments of characteristic polyno- 
mials for random matrix ensembles with unitary symmetry are intimately related 
to semi-classical orthogonal polynomials. By developing the theory of particu- 
lar semi-classical polynomials, it has been possible to characterize these random 
matrix quantities in terms of both discrete and continuous Painleve equations 
[23[2S[TS[3[g[![2g[I3[l3[^ Two methods of devel- 
oping the theory for this purpose have emerged. One has been to use a formulation 
in terms of Lax pairs for isomondromic deformations of linear differential equa- 
tions 23 . The other has proceeded via a theory of ladder operators for orthogonal 
polynomial systems |11L 112] . It is the purpose of this paper to compare these two 
methods as they apply to the particular discontinuous semi-classical weight 

(1.1) w(x) = (1 - C6>(x - t))(x - t^x^e-*, 

with support on M. + , £ < 1, where 8(y) denotes the Heaviside function 9(y) = 1 for 
y > 0, 9(y) — otherwise. 

1.2. A random matrix context. The weight (jl.ip is relevant to the Laguerre 
unitary ensemble LUE^ specified by the eigenvalue probability density function 

1 N 

(1.2) ci[ x ? e ~ xi n (^-^) 2 ^i>o. 

1=1 l<j<k<N 

For a a non-negative integer, this is realized by so called Wishart matrices X'X 
where X is an n X N (n > N, a = n — N) complex matrix of independent standard 
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complex Gaussian matrices (see e.g. [14]). To see how (|1.2j) relates to consider 
the random matrix average 

N 

(1.3) ( - C0(xi - t))\x, - t\»xf 



lue« 

This average is proportional to the multiple integral 

N 

(1.4) Ajy = -r^r / dxf- I dx N T\w(xi) Y] (x k -Xj) 2 , 



A w = — / dxf- I dx N Y\_ w ( x i) H (^fc-^jy 

' ^° ^° i=l \<j<k<N 

where w{x) is given by Introducing the moments 

/>oo 

(1.5) jj, n := I w(x)x n dx 



it is an easy result that 

(1-6) A N = det[^j + fc_2]j-,fc=i,...,iv. 

Moreover, the {A„} can be calculated through a recurrence linking the orthogonal 
polynomials {p n (x)} associated with w(x). 
Introducing the bilinear form 

(1.7) (/,<?) = J f(x)g(x)w(x)dx 

the orthogonal polynomials are specified by the requirements that p n be a polyno- 
mial of degree n and 

(1-8) (pi,Pj) = s i,j- 

By way of application of the Gram-Schmidt process, the weight function and its 
associated support completely specify the coefficients of each polynomial. Fur- 
thermore the orthonormality of the polynomials implies the three term recurrence 
relation [30] 

(1-9) a„+ip„+i = (x- b n )p n - a„p„_i. 

It is the coefficient a n herein which links with A„. Thus one has 

(1.10) a n = . 

Also relevant is the fact that the multiple integral 

(1.11) D N (yi, y 2 )[w(x)\ := — I dxi ■ ■ ■ j dx N \\w{xi){yi - xi)(y 2 - xi) 

l<j<k<N 

can be expressed in terms of the polynomials {p n } according to the Christoffcl- 
Darboux summation (see e.g. |14j ) 

n 10 n n ( u / \i Ajy PN+i{yi)PN{y2) -PN+i(yi)PN(y2) 

(1.12) D N {yi,y 2 )[w{x)\ = . 

7^7^+1 m - 2/2 

In random matrix theory the average (|1.3[) in the case p, = is the generating 
function for the probability that there are exactly k eigenvalues in the interval 
(0,t). In the case £ = it gives the moment of the modulus of the characteristic 
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polynomial Il^i^ — ^O- Moreover, selecting the coefficient of ( in (jl.lj) . taking 
/j = 2 and substituting in shows 

(1.13) D N (s, s)[9(x -t)(x- tfx a e- x ] 

= w d Xl --- dx N Y[( Xl -t) 2 ( Xl - S ) 2 x?e~ x ' l[(x 

After multiplying by 



\Xk Xj j 

1=1 j<k 
-t-sf„ -A2 



this integral can be recognized as being proportional to the joint probability density 
function of the first and second smallest eigenvalues, denoted s and t respectively, 
in the Laguerre unitary ensemble |18j . According to (|1.12p we can calculate (|1.13|) 
in terms of quantities relating to the polynomial system for the weight (|1.1|) . 

1.3. Main results. To present our main results requires some facts from the 
Okamoto r-function theory of the fifth Painleve equation 26j. The fifth Painleve 
equation, 

(1-14) y"=(l + ^ Yy'f-y + ^Uy ■ " 2 



+ 



2y y-lV v " t' i 2 V y 

"32/ . a 4 y(y + 1) 



t y-l ' 

is one of the six non-linear differential equations identified by Painleve and his 
students as being distinct from classical equations and having the special property 
that all movable singularities are poles. The Okamoto r-function theory relating 
to (|1.14[) is based on the Hamiltonian system with Hamiltonian specified by 

(1.15) tH = q(q - 1)V - ((«2 ~ V X ){q - l) 2 - 2{v x + v 2 )q(q - 1) + tq^jp 

+(v 3 - «i)(«4 - vi)(q - 1). 
Here v\ , . . . , v 4 are parameters constrained by 

V\ + V 2 + 1'3 + V4 = 0. 

An essential feature of the theory is that eliminating p in the Hamilton equations 

, _ ,_ OH 

^ dp ' P dq 

(the dashes denote differentiation with respect to t) gives that q satisfies (|1.14|) with 

«i = i(w 3 - w 4 ) 2 , a 2 = -\{V2 - vi) 2 , a 3 = 2vi + 2v 2 - 1, a 4 = -^. 
The first of our main results may now be stated. 
Proposition 1.1. Let v%, . . . ,v 4 be such that 

v 3 - v 4 = -fx, v 3 - vi = n + a, v 3 - v 2 = n, v 2 — V\ = a, 

and so 

a 2 a 2 . , 1 

ax = —, a 2 = — — , a 3 = -{2n + a + 1 + /x), a 4 = -- 
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In terms of the coefficients {a n ,b n } in the three term recurrence i f 1.9\) and the 
parameters a,fi,t of the weight let 

(1.16a) 6„ = b n — 2n — 1 — a — fi — t 

n-l 

(1.16b) «n=(n + f)* + ^-X) 6< - 

We have that the Hamilton equations are satisfied by 
(1.17a) q = — — 

(1.17b) p = 

and so 

(1.18a) 9 n 



t(6 n + t) 



t 



(1.18b) K n = t(n + a + fi/2-pq). 

This characterization is made unique by the specification \S.19\) of the small t ex- 
pansions of 9 n and n n . 

We can also express n and K n in terms of a solution of the fifth Painleve 
equation (|1.14|) with parameters different to those given in Proposition ll.il For this 
we require the fact |26j that (|1.14p is formally unchanged upon the transformations 

(a, [3, 7, S) (->■ (-/3, -a, -7, S) y ^ - 

V 

Corollary 1.2. Suppose v\, ... ,1)4 in H.15\) are such that 
V2 — vi = —pi, V3 — U4 = a, v-2 — V4 = n + a + 1, ua — i>3 = n + 1, 
and furthermore t 1— » —t, so that q satisfies with 

a-i = — , a 2 = — — , a 3 = 2n + a + 1 + \l, on = —~. 

(note that mapping f h> —t in jl-14\ l * s equivalent to mapping a 3 t— > —as). In this 
case the Hamilton equations are satisfied by 

a 

(1.19a) q = 
(1.19b) p = 

and so 

(1.20a) 6 n : 



+ t 

n + t)(n n - fi/2 + 6 n (l + 2n + t + a + fj, + 

t6>„ 



l-q 

„«n * 2 *(*+! -?(2n + a + l + ^)) 

1.20b K„ = ^ + ^ 2 ^ 2 -LL 

(1-q) 2 l-q 
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2. Differential equations for orthogonal polynomial systems 

In the classical theory of orthogonal polynomials, the study of differential equa- 
tions satisfied by orthogonal polynomials has a long and distinguished history [29j . 
Under certain conditions [8], given a system of orthogonal polynomials, the deriva- 
tives of the polynomials may be expressed in terms of a linear combination of at 
most two polynomials of the same system I22L \6\ lllj . To describe these differential 
equations, we parameterize the coefficients of the polynomials according to 

(2.1) Pn{x) = ~f n X n + ln,\X n ' X + ...+ 7„,„. 

In terms of this parameterization, it is clear from (|1.9p that 

7n-l 

a n = 

In 

, _ 7n,l _ 7n+l,l 
In ln+1 

Our starting point is the condition that the logarithmic derivative of the weight 
is rational, and hence we write 

(2.2) S '■"»<*> = f 

where W and V are polynomials in x. This condition implies that the moments 
given by fll.5[) satisfy some recurrence relation. In particular, if we define the 
Stieltjes function by 



f ( X ) = jr» kX -*-i= [ ^l ds 

It x — s 

7. — n d J 



k=0 

then this recurrence is equivalent to / satisfying the holonomic differential equation 

(2.3) W^f = 2Vf + U 

where U is some polynomial in x of degree less than the degree of V . 

We now define the associated polynomials </>„_i and associated functions e n by 
the equation 

(2.4) fp n = e n . 
Explicitly 

(2.5a) e„ = / M£iw s )ds 

JlX-S 

(2.5b) cpn-x = / w{s)ds 

J I s - X 

showing that <f> n -\ is a polynomial of degree n—1 and e„ is meromorphic at x = oo. 
On this latter point, by orthogonality, 

(2.6) c ~ 7„" 1 ^""" 1 

as x tends to oo, while upon multiplying (| 1 .9(1 by / it is also clear that the sequence 
of associated functions {en}^L satisfies (|1.9|) . Using (|1.9p and (|2.1| . we have the 
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large x expansions, 
(2.7a) p n = 

7 » U" - ^ E 5 * + ^ 2 \ E E b * b i E a H + (* n ~ 3 ) 
y j=o y«=i j=o »=i y 

(2.7b) e n = 

(n I n i n+1 

a.—! + x— a x: 6* + *- n ~ 3 E E & ^ + E °? ] + ° (^ 3 ) 
i=0 yi=0j=0 i=l 

Equating (fp n )p n -i with (fp n -i)p n shows <^ n _ip n _i+e„p n _i = <^ n _ 2 Pn+e„. This 
and (|2.7p gives 

(2.8) Pn^n-\ - Pn-l^n = Pn-l<t>n-l ~ Pn<t>n-2 = • 

a n 

The above describes the notation and set formulae to be used in coming sections. 

2.1. Derivation via recurrence relations for moments and isomon- 
odromy. In this section methods are outlined for obtaining differential equations 
satisfied by the orthogonal polynomial system based on the existence of a recur- 
rence for the moments of the weight function, namely (|2.3p . This coupled with 
the theory of isomonodromic deformations allows us to construct the differential 
equations that govern evolution of the polynomials in both x and t. 

Theorem 2.1. The orthogonal polynomials corresponding to a weight w satisfy 
the differential equation 

(2.9) W(x)—p n (x) = (fl n (x) - V(x))p n {x) - a n Q n (x)p n - X {x) 

ax 

where Q n and 0„ are polynomials given by 

(2.10a) 8„ = W {^n-^Pn ~ Pn ~^ en ^J + 2Ve riPn 

(2.10b) fi„ — a n W ^ e "-i^~P" _P " _1 (lL e ™) + a nV( e n-lPn + e„p„-i). 

Proof. First we note from (|2.3p and (|2.4|) that 

w d_ U„-i\ _ 2V0n-i _ u= Wen _ w d_ fen 

da; y p n ) p n p n dx \p 71 

Multiplying through by p n shows (|2.10a|) can be rewritten 

(2.11) 9 n = W [Pn^n-l - tn-l-^Pr^) ~ 2V(p n -i Pn - Up 2 n . 

This tells us that Q n is a polynomial, while (|2.10ap bounds the degree. Explicitly, 
examining the x — > oo behaviour, namely (12. ip and (|2.6|) . shows 

(2.12) deg6„ < max(degV7-2,degV-2,0) 
as noted in 
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Using (HU) and (|2"TT|) we find 

a n (Pn-l<Pn-l ~ Pn<t>n-l)®n = 

W yp n ^(/) n ^i - ^-l^Pr^j _ 2V(f) n -iPn - Up n . 

By appropriately grouping terms divisible by <f> n -i and p n on opposite sides, we 
define the polynomial Q n to be the common factor according to 

(2.13) Pn<t>n-l^n = 0n-l ( 0>n®nPn-l + W ~j~P™ + ^™ 

= Pn ^a„6„</> n _2 + _ V0n-1 - E/Pn 

The first expression in (j2~T31) is equivalent to (J3HJ) provided (|2.10b[) can be verified. 
For this purpose we use (|2.10a[) in (|2.13|) to obtain 

Pn Pn Pn Pn 

By rearranging (|2.8|) . we find a n e n p n -i — a n p n e n -i — 1, which we use to remove 
occurrences of p n -i, giving (|2. 10b)) as required. Examining the x — > oo behaviour, 
by using JUT]) and JH]) in (|2.10b|) . shows 

(2.14) degfi n <max(W-l,y) 

which again appears in |23j . □ 



The origin of this theorem can be traced back to the work of Laguerre 
and has since been revisted by contemporaries [H [3 I23j . The theorem provides 
a mechanical way of determining the differential equation satisfied by polynomials 
provided one knows the rational logarithmic derivative. One need only expand 
(|2.10a[) and (|2.10b[) to polynomial orders using (|2.7b[) to produce a parametrization 
of the differential equation satisfied by the polynomials in terms of the a„'s and 
b n 's. A simple application of (|1.9[) gives us an expression for the derivative of p n —i, 
which is but one column solution to a 2 x 2 linear differential equation in x. The 
following corollary provides us with another solution. 

Corollary 2.2. The function e n / w satisfies (|2.9|) . 

PROOF. Consider the derivative of fp n in terms of e„ and <p n . According to 



W-^-fp n = W^-4>n-i + W^-e n . 
ax ax ax 

On the other hand, use of (|2.9[) and (|2.3[) shows 

W^-f Pn = f[W^-pn)+ Pn W^-f 

ax \ ax J ax 

= (fl n - V)fp n - a n O n fp n -i + 2Vfp n + Up n 
= (Q n + V)(j) n -i - Cl n Q n (j) n -i + Up n 

+(f2„ + V)e n - a„6 n e n _i, 

where in obtaining the final equality (|2.4[) has also been used. By cancelling out the 
derivative of 4> n -i as calculated from the first expression for Cl n from the previous 
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proof, we deduce that the derivative of e„ is given by 
Hence 

w d_en _ wW£e n -e n W±w 



dx w w 2 

(On + V)e n ~ QnQ«e«-i - 2Ue» 

where use has also been made of (|2.2[) , as required. □ 

As a linear system, we have two linearly independent solutions. As mentioned 
above, both {p n } and {e n /w} satisfy (|1.9|) and (|2.9p . telling us that p n -\ and e n _i 
satisfy 

W—y n = a n e„_iy„ + (O n _i - V - (x - 6„))y„_i. 
ax 

Hence, the matrix 



(2.15) Y n 



Pn w 

Pn-l . ^T" 



satisfies the matrix differential equation 

(2.16) -^y„ = <y„ 

where 

> _ J_ / - v -a n e„ 
" W" \-a„0„-i n„_i - F - (x - 6 n _i)6 n _i 
Because and {e„} satisfy (|1.9p . y„ also satisfies 

(2.17) Y n+1 = M n Y n 
where 

c — b n a. 



\ i o 

Lemma 2.3. The polynomials 0„ and tt n satisfy the recurrence relations 

(2.18a) W + al +1 Q n+1 - a*0 n -i = (x - 6„)(0„+i - fi„) 

(2.18b) (x - 6 n _r)9 n _i - (x - &„)©„ = fi„ - Q n+ i 

Proof. The equation (|2.16p and (|2.17|) gives two ways of calculating -^Y n+ i. 
The consistency may be written as 

dM n 

M n An ~ An+lM n + — — = 0. 

dx 

This is an identity on the bottom two rows, however, in the first row, the consistency 
relation yeilds (|2.18aj) and (|2.18b|) □ 

This lemma can be found in Magnus [23 . Using (|2.8p . we have 

detF„ = . 

a n w 

In general, for an equation of the form (|2.16p , we see from (|2.15p that 

— det Y n = Ti£/ n det Y n 
dx 
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and so 

-r— = ^TT7 = Tr^n det y„ 

aia n w wa n W 

giving the additional relation 

(2.19) (x - b n )e n = n n+1 + n n . 

This also implies (|2.18b[) . It further gives us a new parameterization of &/ n , given 

by 

. _ J_ /f2„ - V -a n @ n 
W \a n Q n -i -Q n — Vj 

as first derived in [23] . 

Another useful relation comes from the multiplication of (I2.18a|) and (|2.19|) . 
which gives 

(2.20) we n + a 2 n+1 e n+1 e n - a 2 n Q n <d n -i = n 2 n+1 - nl . 

Summing over n, given that f2o — V, shows 

n-i 

(2.21) n* -alQ n e n ^ = v 2 + wJ2 Q i- 

8=0 

The roots of W are now the poles of sf n . If {xj} is the set of poles of s$i, then 
we may write (|2 . 16[) as 

(2.22) ±Y n = , <Yn = [ J] _£%. ] r„ 



. _ 1 (it) - V(xi) -a n Q n (xi) 



where 

1 !"'(..-,) V flnQn-lfe) -fin(l«) - Vfo)/ 

This now places the differential equation into the context of isomodromic deforma- 
tions. In general, any solution to (I2.22p is going to be multivalued, with branch 
points at {re,;} (one possibly being oo). Hence, by integrating around a path, say 
p : [0,1] — > C \ {xi} where p(0) — p(l), the multi-valuedness can be expressed 
through the equation 

Y(p(0)) = Y{p{l))JZ p 

where is referred to as a monodromy matrix. The set of monodromy matrices, 
{^ p }, forms a representation of the fundemental group, 7i"i(C \ {x,}). This is 
depicted in figure [1] 

By construction, one solution of (|2.22[) involves the polynomials, which are en- 
tire, and the associated functions. The goal of monodromy preserving deformations 
is to describe a family of linear problems of the form (|2.22[) that share the same 
representation. A natural choice of deformation parameter turns out to be the poles 
of An, giving rise to the classical result known as the Schlesinger equations, given 
by 

a* _ Ia^M 

OCtj (Xi — Ctj 

(2.23b) p = -r&lM. 

3^1- 




Figure 1 . The fundamental group of the complement of a set of 
three poles in C. 



We shall assume that just one of the poles depends on a variable, t, which shall 
become the deformation parameter. In the case that we have just one parameter 
that needs to be deformed, we have that there is a matrix 

BY 

(2-24) m n {x,t) = ^Y-^ 

The form of this matrix, as implied by (|2.23p , is given by [13] 

(2-25) SS n — &oo,n ~ ^ : TTT ■ 

X X ^ Cjl 

i 

By examining the large x behaviour of p n and we deduce that 3$oc n is given 

by 

/_L2Za o \ 

(2-26) ^oo,„= 7 " 91 x a7 „_ x • 



\ Tn-l dt / 

This gives us two linear differential equations for the one system, which must be 
consistent. Hence we have the condition 

— —Y n = — — Y n , 

dt dx n dx dt ' 

which is equivalent to 

d d 

(2.27) sg n m n - 3§ n srf n + —y n - — SB n = 0. 

at ox 

This completely determines the differential equation for the orthogonal polynomials 
and associated functions in t. 

2.2. Ladder Operators. An alternative approach, developed by Chen and 
collaborators [111 I12L l4j. is that of the ladder operators. We shall assume that the 
weight satisfies the same logarithmic differential equation (|2.2p . and hence that the 
corresponding moments satisfy (|2.3p . However, this approach typically concerns 
the monic versions of the orthgonal polynomials, given by 

P n = —Pn- 

In 
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In order to make comparisons to the previous section, we will deal primarily with 
p n rather than P n . Now p n , being a polynomial of degree n, when differentiated 
can be expressed as a linear combination of {pj\j=o,...,n—i, 

d n— i 

(2.28) — p n = a n ,iPi- 

x i=0 

We may reduce this, via the use of (|1.9[) . to a differential equation specified by the 
following theorem of Bonan and Clark [7] and Bauldry [5] . 

Theorem 2.4. The orthogonal polynomial system {p n }, defined by the weight 
w, satisfies 

(2.29) 4~Pn = SnPn + a n A n p n _ 1 

ax 

where 

V ' Ji\x-vJ\ W(x) W{y) ) yy > 



(2.30b) B n 



f f p n (y) Pn -i(y) \ (2V(x) 2V(y)\ 

an l { x-y ) [wm ~ urn) w{y)dy - 



PROOF. Beginning with l|2.28[) . we may use orthogonality and integration by 
parts to find 

a n ,z = J^p'„(y)Pt(y)'w(y)dy = - J^Pn(y) (p'i(y) ~ -^jj^— J w{y)dy. 

However, since p' t is a polynomial whose degree is less than n, this term must be 
destroyed by orthogonality, leaving 

f 2V 

= ~ J r -^Pn(y)Pk(y)w(y)dy. 
The derivative of p n can therefore be written 

/ n ~ 2V^(y) 
^2Pi(x)Pi(y) w ^ w(y)dy. 

However, if we replace 2V(y)/W(y) with 2V(x)/W(x), this would vanish by or- 
thogonality, hence we may add it to obtain 

p'M = J My) (|>(*)ft(y)) (fg - Hy)*y. 

By exploiting the Christoffel-Darboux summation (|1.12|) and pulling out the poly- 
nomials in i, we arrive at 

Pn (x) = a n p n . 1 {x)J i (— — j (_ - W (y)dy 
f f Pn (y) Pn -i(y) \ ( 2V(x) 2V(y)\ 

" a " p " (x) X { x -y ) [wixj ~ my)) w{y)dy > 

which is 1(235)1 . □ 
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This allows us to define the ladder operator, 
which has the effect 

L n ,lPn = AnVn-l- 

Lemma 2.5. The terms A n and B n satisfy the recurrence relations 

2V 

(2.31a) B n+1 +B„ = ( x -b n )A n - — 

(2.31b) {B n+1 -B n )(x-b n ) = a 2 n+1 A n+1 - alA n ^ + 1 

Proof. Using (|2.30b|) . 
o , o f ( Pn(y){a n Pn-i(y) + a n +iPn+i(y)) \ ( W{x) 2V(y) \ 

+ "I, { —« ) [wM - wW)) aiy)i " 



= [x — b n )A n — -. 
By consistency of f|2 .29[) with (| 1 . 9|) we obtain the second required expression. □ 

By multiplying and rearranging (|2.31a|) and (|2.31b[) we obtain 
2V 

(2.32) B n+1 - B n - — (B n+ i - B n ) = a n+1 A n+ iA n - a n A n _ x A n + A n . 
Hence, by summing over n and appropriately evaluate initial conditions, we obtain 

2V n ~ 1 

(2.33) B 2 n - —B n - alKK-x = - E A - 



Note the structual correspondence of (|2~32"|) and (|2T55|) with (f2T2TJ|) and (|2T2Tj) re- 
spectively. 



3. Derivations of Py 

We now turn to the polynomials specified by (|1.7p with the weight specified by 
(II. ip . The above formulae for the derivatives in x should be considered as partial 
derivatives. The formula (|1.5[) for the moments is a hypergeometric integral, which 
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may be evaluated to give 



(3.1) 



/'■A' 



(C- l)T(l + k + a + n)iFi 



—k — a — fi 



-t + 



(C - 1) sin(7r(a)) \ T(k + 1)I> + 1)! 



t 



sin(7r(a + fi)) J T(2 + k + a + fi) 

l+k+a+fj, F [ I + k + a 

2 + k + a + /i 



= Ci(C,/i, fc, a)iFi 



- t 



-fx 



where 



— k — a — /! 



d((,n,k,a) = (l-C)r(l + fc + a + ^) 

(C-l)sin(Tra ) ^ T(fc + l)T(p + 1) 



C 2 {C,n,k,a) 



(-iy 



sin(7r(a + /i)) J T(2 + k + a + fi) 



and iFi is the confluent hypergeometric function. We seek the corresponding dif- 
ferential equations satisfied by the orthogonal polynomial system, as implied by the 
theory of JO] and 32~2l 

To derive the differential equation satisfied by the orthogonal polynomials from 
the theory of ij2.ll we remark that since the factor of (1 + C,0(x — t)) plays the 
role of a multiplicative constant almost everywhere, the logarithmic derivative of 
w coincides with the logarithmic derivative of (x — t) tJ/ x a e~ x almost everywhere. 
Hence we write 

x(x — t)— — = {—x 2 + (a + /i + t)x — /it) 

w 

where = is to be interpreted as equals almost everywhere, and so independent of £ 

W = x(x-t) 

2V = -x 2 + (a + fi - t)x + /.it. 

Recall that the form of the logarithmic derivative is the essential ingredient in both 
theorem 12. II and theorem [ 



3.1. Recurrence of moments approach. Now that W and V have been 
defined, determining the differential equation satisfied for this particular family of 
orthogonal polynomials is simply a matter of applying theorem 12.11 

COROLLARY 3.1. The polynomials p n corresponding to the weight (|1.1[) satisfy 
the differential equation 

(3.2) dJ Pn )=\^ + ^ + ^\( n Pn 

\Pn-l I 1 X X~t\ \Pn-l 
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whe 



1 /k- _ M _ n ft 

= ~ a _ , _ ^ 

nf = - ( ( n + 2) 1 ~ Kn a n (O n + t) 

t\-a n {6 n -i+t) K n -(n + a+f)t 

/; '!! 1 



2n — 1 — a — (i 



(»+f) 



i 1 2 , 7n,l 

t + a n -{ '-. 



In 

Proof. By way of application of (|2.10b[) and (|2.10ap using (|2.7|) one obtains 
for the explicit form of O n and n , 

x 2 ( 2n + a + ji + t\ 2n + [it 2 7n,l 



= ~T + H 2 J 2 7 „ 

9„ = — x + 2n + 1 + a + fx + 1 — b n , 

which we subsequently decompose into the form seen above. □ 

We also require these equations to be written in terms of the K n and 6 n alone. 
For this, we note that recurrence relations for 9 n and K n are implied by (|2.18ap and 

snm . 

COROLLARY 3.2. The associated functions, 9 n and K n , satisfy the recurrences 
(3.3a) K n+ i + K n = -0 n {On + t + 2n + a + 1 + /i) 

ft ft K 2 - 

u n u n—l "fn 4 



(3.3b) 



?n + t Q n -i +t (n n - (n + a + f ) t) («„ - (n + § ) t) ' 



Proof. The relation (|2.18b|) is equivalent to (|3.3ap when one uses the defini- 
tions of n and f2 n in terms of 9 n and K n . Evaluating (|2.21[) at x — and x = t 
shows 



(3.4) <£Mn-l = « 2 



^2 

4 



(3.5) fl »(t + * B )(t + *»_ 1 ) = (^- (2 " + 2 ; + /i)t ) (/C-^l^* 

respectively. The ratio of these identities is (|3.3b[) . □ 



The relation (|3~3b"l) may be used to eliminate the occurence of in 
Note that (|3.2p is the form of (|2.22p . where only one of the poles depends on t. 
Hence, the evolution in t is governed by (|2.24[) . In this regard, the derivations 
of the time derivatives for 6 n and K n and the methods of Forrester and Witte 
[18) contrast with the methods of Basor and Chen [4]. Once the derivatives in x 
are found, one may apply the theory of isomondromic deformations [21] to obtain 
appropriate derivatives in t. This approach, as seen in [18] extends the evolution 
of the orthogonal polynomials in the t direction via the following result. 
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COROLLARY 3.3. In addition to S3. 2\) . the othogonal polynomials satisfy 
(3.6) di < Pn\_f^ *t\{Pn 

where 



Pn-l J X - t I \Pn-l 



m _±(e n + t o 

2t \ -9 n ^ - t 

Proof. This almost directly follows from the corollary of the Schlesinger equa- 
tions, (|3.2[) and ()2 . 25[) . We use the fact that in the context of orthogonal polyno- 
mials has the explicit form (|2.26[) . By equating the residues of the left and right 
hand side of the compatibility relation (|2.27p at x — oo, the diagonal entries reveal 

(3.7a) = l+ °n 

In t 

(3.7b) = 1 + ^1 

7n-l t 

while the off diagonal entries are 0. This gives the required form for 2$ above. □ 

One may easily calculate the derivatives in t of 9 n and n n via the compatibility 
of (|2~T1) and (l2~m 

(3.8) d t s# - d x SS + - = 0, 

to define the evolution of 9 n and K n . Using the above recursion relations allows 
one to express the derivatives of 9 n and n n in terms of themselves. Alternatively, 
using the general framework of |23j , the derivatives of a n and b n are expressible in 
terms of the functions n and f2 n evaluated at the movable finite singular points 
of (|2.16p via the expression 

(3 ' 9a) d7 lna " = 2^ wW) Jt Xr 

(3 ' 9b) dt lnbn = 22 WW) dt Xr > 

r—1 v ' 

where in the case of (|2.16[) . m = 1, and the only point is x\ = t. This leads to the 
equations 

(3.10a) = 2 + fo„_ 1 -6„ 

a n at 



(3.10b) t-ji = 



. db„ 
df 

We know (|3.9|) is equivalent to (|3. 10[) . Using (|3.8[) and (|3. 10[) in conjunction with 
(|3.3[) to eliminate occurences of o? n and n _i gives a differential system for {#„, n n }. 

COROLLARY 3.4. The associated functions, 9 n and K n satisfy the coupled dif- 
ferential equations in t 

(3.11a) t-§- t 6 n = 2k„ + {2n + a + l + fi + t + 9 n )9 n 

0k„ _ ( 1 | 1 \ „2 



(3.11b) t^ = ^ g _L_ + ^_j^ + ^2n+a + /x + l-(2n + a + M) 5 ^ 
- (n 2 + (n + 1 ) (a + //)) t - + ( n + H) ( n + Q + g) ^ 
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Proof. This simply follows from the evaluation of (12.27| . The first relation 
follows from (|3.7|) . namely 



— t>n-l — Vn- 

The two other relations that arise are 

t6' n = 2n n + (2n + a + l + f i + t + 8 n )9 n 
tn n = K n — a n {Q n 6 n — i). 

The first of these is (|3.11ap . By using (|3.4[) to eliminate a n and (|3.3b[) to eliminate 
6 n -x, one obtains (|3.11b[) □ 

Now that one has the derivatives of K n and 9 n , the remaining task is to find 
the transformation which allows them to be identified as the Hamilton equations 
for a Painleve V system. But before doing this, we want to show how differential 
equations equivalent to the coupled system (|3.1ip can be derived from the formalism 
of §2.2. 

3.2. Ladder operator approach. We want to specialize (|2.29[) to the weight 

m- 

PROPOSITION 3.5. With the weight the coefficients A n and B n in (|2.29j) 

are given by 

. _ Rn 1 — Rn 

A n — 



B„ 



x — t x 

r n _ n + r n 
x — t x 



where for a > 1 



Rn = a I {t _ y) dy 



, w{y)Pn{y)Pn-l , 

al (t-y) ^ 



Proof. We note that 

2V(x) 2V(y) fa fx 



W{x) W{y) \ (t-x)(t-y) xy 
and hence the integrals that define A n and B n simplify to 



(x - y), 







- 1 " ^ 1 Pniy)2 '- (t- X )(t-y) -w) lr{!M!l 



Bn = / Pn(y)Pn-l(y) ( -77 777 7 ~ — ) w(y)dy, 

V \t-x)(t-y) xy; 
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or equivalcntly 



An = — I (l-Ce(y-t))(y-tr-Ve- y p n (y) 2 dy 
x — t ' 







X 



± I (l-C0{y-t))(y-t) a y> t - 1 e-Vp n (y) 2 dy 



(i 



B n = —J {l-Ce{y-t)){y-tr-\fe-y Pn (y)p n ^{y)dy 
x — t ' 





1 r°° 

— (1 - C% - t))(y - try^e-^p^pn-^dy. 
x Jo 

Now define 

poo 

Rn = a {\~^{y-t)){y-tr-\fe-y Pn {yfdy 



[°° w{y)p n {y)\ 
a / — 77 \ — dy 



o (t-y) 

POO 

r n = a (l-Ce(y-t))(y-t) a - 1 y^p n (y)p n . 1 (y)dy 



o 



n 



/■°° w(y)p„(y)p„-i 



We apply integration by parts, orthogonality and the known value of w(x, t) at 0, t 
and oo, to express the second part of the integrals in A n and B n in terms of R n 
and r n respectively, giving 

a, = - i r r^fiMM _ ^.tKo,)) dy 

z-f a:,/o V t-y J 

o r n n f°° ( aw(y,t)p n (y)p n ^ 1 ,dp n \ 

= ^-x7 Q v + ^,iK-i(y)^Jd y . 

Using orthogonality and the expression 

— ^ = — p n -i(y) + lower order terms 

a n 

gives the stated formulas. □ 

Now that the form of A n and B n are known, the differential equation satisfied 
by the polynomials can be written 

d 



da; \Pn-i 



n+r _ r„ ( l—Rn i Rn 

Pn \ [ a x-t Un \ x ~T x-t J \ ( Pn 

( l-R n -t Rn-i \ r„ n+r n 2V I [p„_i 

in \x x-t ) x-t x W / y 



where the second row is a consequence of (|1.9|) and (|2.31al) . 

Lemma 3.6. The functions R n and r n satisfy the recurrences 
(3.12a) r n +i+r n — a — R n {p + a + 2n + 1 + tR n — t) 

(3 12b) RnRn—i _ r n (r n Cl) 



(i^-l)^-!-!) (r„ +n)(r n + n + p) 
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Proof. The residue in x of (| 2 . 3 1 a|) at t and oo using these definitions for A n 
and B n in terms of R n and r n shows 

(3.13) b n = 2n + 1 + a + fi + tR n 
r n +i + r„ - a = R n (t - b n ) 

which gives ()3.12ap . The evaluation of the result of multiplying (12.31bj) by x 2 (x — f) 2 
at and t reveals 

(3.14) r„(r n -a) = a 2 n R n -iR n 

(n + r n )(n + fi + r n ) = a 2 n (R n - l)(R n -i - 1) 
giving (|3.12b[) in an analogous manner to the previous section. □ 

In addition to (|3.13[) . there is a further relation obtained by eliminating R n -i 
from ()3.14j) by using (I3.12bj) . giving 

/ 01c x 2 (r n -a)r n (n + r n )(n + fi + r„) 
(3-15) a n = . 

K n tin — t 

Lemma 3.7. The recursion coefficients of (jl.9p satisfy the differential equations 

(3.16a) = R n - 1 -R n 

a n at 

(3.16b) — j- = r n - r n+ i 

Proof. We take the derivative of (|1.8p in the case i = j = n to see that 

d f 00 

= -j£ J p n (y) 2 w(y)dy 

with respect to t. We recall that a is a non-negative integer by assumption. Hence, 
w(x, t) is continuous at x — t and w(t, t) = 0, and so 



= ^ ( / <ly ■ / <!'/ j ic{y)i>„\y) 2 



lim — lim I w(x)p n (x) - 

x — >t~ x — >t+ J 

dy+ / 2p n —-w(y)dy 



f dy - J dy^j ^ (p n (y) 2 w(y)) 



y — t Jo dt 

where we have used 

dw(x,t) aw(x,t) 

dt x — t 



Using 

gives us that 



d Pn In n , i j 

-—— = — jn X + lower order terms 
dt 7„ 



Rn = 2^ 

In 



and hence 



2 da n _ Q ( 7 n _i \ 1 d ( ln -i \ _ < £h-i _ <£L =R R 



a n dy V 7n / dt \ j n J 7„_ x 7., 
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Similarly, differentiating (|1.8|) in the case of i = j + 1 = n with respect to t 
shows 

/>oo 

{i' n y n + in,iy n ~ l )pn-i{y)&y 





= — I" J n y n Pn-My)dy + ^ i' n , 1 v n - 1 p n -Mv)dy. 

In Jo In Jo 

We use the fact that 

7„a;™ = p n (x) — 7n,i#"~ + lower order terms 
in this expression to obtain 

oo / 7 



/ / 7n,l7n\ n-1 - 

(7n,i - Jy Pn-iwdy 



In ln,\ln ~ 7n7n,l 



since 



7n-l Jo \ In 



7n,l 



^ _ In ,1 _ 7»+l,l 



7n 7n+l 

the equation (|3.16bp follows. □ 

Theorem 3.8. The coefficients, R n and r n , satisfy the system of differential 
equations 

(3.17a) tR' n = 2r n -a + R n (tR n + 2n + a + fj,-t) 

(3.17b) tr' n = ( ^- 2R " \rl-n{n + a ) U " 



Rn(l — Rn) J ™ 1 — Rn 

+(2n + a + fi)r n H _ — 

PROOF. To obtain the first equation, note that 

^ = i?„ + ii^ = r n - r n+1 

which gives (|3.17a[) under the substitution of r n+ i in accordance with (|3.12a|) . 

To obtain (|3.17b|) differentiate (|3 . 1 5|) . Knowing a' n in terms of a n , r n and R n 
from (|3.16ap removes a' n while we may use (|3.14[) to remove the remaining instances 
of a n . Using (|3.17ap and ()3. 12b[) are used to eliminate R' n and R n -i to obtain an 
equivalent reformulation of (|3.17b[) . □ 

We observe that the differential equations of corollary 13. 21 and theorem 13.81 are 
identical upon setting 

t + n 



t 

K n ( (1 

— In 



t V ' 2, 

This demonstrates an equivalence between the characterization of the polynomial 
system corresponding to (jl.lj) as implied by the method of isomonodromic defor- 
mation, and the method of ladder operators. 
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3.3. Main results. Recall that our remaining task is to relate the differential 
equations of corollary to the Hamilton equations for a Painleve V system. The work 
in [18j and [4] both provide clues regarding relevant transformations. Explicitly, 
they suggest two Mobius transforms 



t + 9„ 

n +t 



both of which send oo to 1, and send and — t to and oo in different ways. 

PROOF of proposition 11.11 We are required to show that q satisfies (I1.14[) 
using (|3.1ip . We first find the derivatives of q in terms of 9 n and «„, 

n +t 



q 
q 



In - t9' n 

9 2 

2n n + 9 n {9 n + 2n + t + a + /j,) 
ff 2 

w it 

9' n ((2n + t + a + (i)6 n + 4/t„) 1 2n n + 9 n 



9 3 9 9 2 

However, using (|3.11a[) we have 

2k„ = t9' n - 9 n {9 n + 1 + 2n + t + a + fi) 

giving 

„ = t(2t + 29 n )9' 2 (2t + 9 n )9' n 1 9\__ 

q 29 3 n (t + 9 n ) 9 2 n {t + 9 n ) 29 n {t + 9 n ) t(t + 9 n ) 

y 2 (t + 9 n ) 2 -9 2 n (t 2 + a 2 ) 1 / t t + 9 n 

+ 29jitT9~) +2 v ^~ 2(1 + 2?l + a + /i) ^T~ 5 

Inverting the expression for q in terms of 9 n gives us 

9 n — r 

q - 1 

q-tq'-l 
{q-lf ' 

and using these expressions show 



ji 



q = i—^ + tt. q - - 



i i \ ,i q 1 , (q-i) 2 (tfq o? 



This is ([01]) where 



q-\ 2q) 1 t t 2 \ 2 2q 

(l + 2n + a + Li)q _ q{\ + q) 
t 2( g -l)' 



IL 2 a 2 
«i = - a 2 = -— 



a 3 = — (2n + 1 + a + [i) a^ = —-. 
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To obtain the corresponding p variable, we remark that the equation for q' , as 
specified by the Hamiltonian in (|1.15p . is linear in p, hence determines p uniquely 
in terms of 9 n and K n . □ 

For the differential equations to uniquely characterize 9 n ,K n , boundary values 
must be specified. For this purpose, we note from (|3.ip that the small t leading 
order asymptotics for p k are 

= c ( 1 tj t mO* - 1) t 2 

1 \ k + a + p 2(k + a + fi)(k + a + a* - 1) 

t3 

6(fc + a + p)(k + a + p - l)(fc + a + p - 2) 
+C2t i+k +a +„ ( 1 _ 1 + k + a t + (1 + k + tt )(2 + k + a) 



2+ k+a+p 2(2 + k + a + /x)(3 + k + a + y) 
(l + k + a)(2 + k + a)(3 + k + a) 3 



6(2 + k + a + ^) (3 + k + a + pi){A + k + a + p) 



From this the determinant (|1.6p that defines A„ may be evaluated to leading orders 
by using the identity (see e.g. [25]) 

n-l 

(3.18) det(T(z k +j)) J . k=Q .._ n _ 1 = Y[T(z k ) JJ (z k -Zj). 

k=0 Q<j<k<n 

In particular, by letting Zk = l+ a + p + k we have 

n— 1 n— 1 

A„(0) = (1 - 0" J] fc! 1] r(l + a + a* + k). 

fe=l fc=0 

Recalling (|1.10p then gives 

a n (0) = n(n + a + p), 

and knowing this (|3.10[) at t — implies 

6„(0) = 2n + a + 1 + 

To determine the rest of the expansion of (|1.6p , we first make use of (|3.18p to 
compute the leading form of the analytic and non-analytic components as 

A n (t) = A„(0) (l + — y — t + 0(t 2 ) + Xn t 1+a+ »(l + 0(t) + 0(t 1+a +n)) 
V a + p / 

where 

Xn = (i - c)"- 1 ((-iy + (c ~ l) sin ^ a) ^ r(M + 1)r(a + ^ + n + l ) 



sin(7r(a + p)) ) (n - l)T(a + p + l)T(a + p + 2) 2 ' 

It follows from this and (jl.lOp that the expansion of o? n is 

? , . np 2 (a + p + n) i7 

a 2 n = n(n + a + p)- ( + ^ 2 1 + 0{t ) 

,,1+q+M r(q + At + n+ i)r(i + ^) ,,i +M+al 
n!r(a + /x)r(a + M + l)r(a + a* + 2) + 1 J ' 
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and subsequently consistency with regards to ()3.10|) demands 

a au(2n + a + u + 1) 9 , 
b n = 2n + a + ^ + 1 1 + - PV 1W , 9 . ; -r t + 0(t 3 ) 

T(q + ^ + n + + 1) 

r(a + /z + l) 3 

Substitution into (| 1 . 1 6[) then shows 

n m \j ■ , a/i(l + 2n + a + ^ 2 3 

y n(*J : — 1+7 — : 7T7 — : — — —. rri + 0{t ) 

a + fi (a + [i - ±){a + /u) 2 (a + 1) 

t i+q+M r(a + m + n + i)r(^ + 1) . 2+Al+a! . 

r(q + /i + l) 3 1 1 

= ,(2n + a + ,) t _ 2nMn + * + ») 
W 2 + a + fi ((a + /i) 2 - l)(a + ^i) 2 v ; 

(a + M + " ~ IjTQx + l)r(a + g + n + 1) 
T(a + [i+ l) 3 (a + n) 

(3.19) +0(i 2+ ^ +a ) 

In particular, the small t asymptotics for q is therefore 

q_6 n + t _ ct a(2n + a + fi + 1) 1 + 0(1?) 
9 n fx fx(a + fi + l)(a + fi — 1) 

+ ^ n»)r(» + a + M+l) +0(f i +M+a)- 
//(/i + a)T(a + [ly 
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